Exact results for a tunnel-coupled pair of trapped Bose-Einstein condensates 
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A model describing coherent quantum tunneling between two trapped Bose-Einstein condensates 
is shown to admit an exact solution. The spectrum is obtained by the algebraic Bethe ansatz. An 
asymptotic analysis of the Bethe ansatz equations leads us to explicit expressions for the energies 
of the ground and first excited states in the limit of weak tunneling and all energies for strong 
tunneling. The results are used to extract the asymptotic limits of the quantum fluctuations of the 
boson number difference between the two Bose-Einstein condensates and to characterize the degree 
of coherence in the system. 
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Macroscopic quantum tunneling is one of the most fas- 
cinating phenomena in condensed matter physics jl],|]]. 
Experimental realization of Bose-Einstein condensates 
(BEC's) in dilute atomic alkali gases has stimulated a 
diverse range of theoretical and experimental research 
activity A particularly exciting possibility is that 

a pair of BEC's (such as a BEC trapped in a double- 
well potential) may provide a model tunable system in 
which to observe macroscopic quantum tunneling. In 
this Letter we show that the canonical Hamiltonian for 
a pair of tunnel-coupled BEC's || has an exact solution. 
The model is also realizable in Josephson-coupled super- 
conducting metallic nano-particles . This connects the 
model to the rich and powerful mathematics associated 
with the algebraic Bethe ansatz || and provides a means 
to study the degradation of the coherence between the 
two BEC's which occurs as the tunnel coupling increases 
(or the size of the BEC decreases) . The exactness of our 
approach means that it is not necessary to resort to ap- 
proximations such as Gross-Pitaevskii mean-field theory 
and the phase-number formulation which become ques- 
tionable in the coherent-incoherent crossover region. 

The canonical Hamiltonian takes the form || 
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— (a[a 2 + aJ,ai). 
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where a\ , a\ denote the single-particle creation operators 
in the two wells and N\ = a\a±,N2 = a\a,2 are the 
corresponding boson number operators. The total bo- 
son number N\ + N2 is conserved and set to the fixed 
value of N. The physical meaning of the coupling pa- 
rameters for different realizable systems may be found in 
Rcf. M. It is useful to divide the parameter space into 
three regimes; viz. Rabi {K/£j « N^ 1 ), Josephson 
(A^ 1 << K/£.j « N) and Fock (N « K/£j). In 
the Josephson region one expects coherent superposition 
of the two BEC's (Schrodinger cat states) to be possi- 
ble whereas in the Fock region the two BEC's will be (in 



some sense) localised in the two separate wells. There is 
a correspondence between ([!]) and the motion of a pendu- 
lum ||. In the Rabi and Josephson regimes this motion 
is semiclassical, unlike the case of the Fock regime. For 
both the Fock and Josephson regimes the analogy cor- 
responds to a pendulum with fixed length, while in the 
Rabi regime the length varies. An important problem is 
to study the behaviour in the crossover regimes. 

The exact solvability of ([j]) which we discuss here 
follows from the fact that it is mathematically equiva- 
lent to the discrete self-trapping dimer model studied by 
Enol'skii et al. jn| , who solved the model through the al- 
gebraic Bethe ansatz. The Lax operator appearing in the 
Yang-Baxter equation which underlies the exact solution 
first appeared in Q. The consequences of this result 
are legion. The exact solution gives a direct method, 
through numerical analysis, to investigate the nature of 
the system in the crossover regions between the Fock, 
Rabi and Josephson regimes. It also provides a means to 
test the limits of applicability of previous approximate 
treatments [p~2|— . Solvability through the Yang-Baxter 
equation raises questions about the probability distribu- 
tion of the energy level spacings pj| . Moreover, the for- 
mulation of the model through the algebraic version of 
the Bethe ansatz opens possibilities for the calculation of 
time-dependent form factors and correlation functions, 
as achieved in [^o| for A/i = 0. 

In this Letter we will show that the exact solution al- 
lows for the analysis of the asymptotic behaviour of the 
energy spectrum in the limits of strong and weak tun- 
neling. These results are used to extract the asymptotic 
limits of the quantum fluctuations of the relative par- 
ticle number between the two BEC's and the degree of 
coherence in the system for the ground state. Our ap- 
proach also makes it feasible to undertake an asymptotic 
analysis of the system at finite temperature which gives 
insight into the relative influence of thermal and quan- 
tum fluctuations (cf. pi)). We mention that these find- 
ings provide a useful tool in the numerical evaluation of 
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the exact solution in the strong and weak tunneling re- 
gions. In contrast to most Bethe ansatz solvable models 
P], the model discused here, as in the case of the Azbel- 
Hofstadter problem , is a quantum mechanical model 
rather than a quantum field theory; i.e., it has a fixed 
and finite number of degrees of freedom. 

Bethe ansatz solution. Following the standard proce- 
dure of the algebraic Bethe ansatz || , we can derive the 
Bethe ansatz equations (BAE). For N total bosons the 
BAE for (1) can be written in the form 



N 



P^a. 
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The parameters 77, oj which naturally arise in the Bethe 
ansatz solution, along with a scaling factor k, are related 
to the coupling constants of ([[J) through the identification 

K = 2nn 2 , Afi — —2Knuj, Ej = 2k. 

Each set of numbers {t) a }„_j which is a solution of the 
BAE defines an eigenstate of the Hamiltonian. For de- 
tails of the derivation of (2) we refer to [ ^OpC| ] . For such 
a solution, the eigenvector has the form 

N 

*(yi,...,v N ) = l[C(v a )\0} 

a=l 

with C(u) — (u — oj + riN 2 )a\ + rj -1 ^ and |0) the Fock 
vacuum. 

For each N we expect N + 1 independent solutions 
of the BAE. Note that in the derivation of the BAE it 
is assumed that v a are distinct for different a. This is 
a result of the Pauli principle for Bethe ansatz solvable 
models as proved by Izergin and Korepin p3| for the 
one-dimensional <5-function interaction boson gas. Their 
result may be adapted to the present case, which plays 
an important role in our asymptotic analysis of the Bethe 
ansatz solutions. From the BAE, we may derive the use- 
ful identity 
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a=l a=l p=m+l M 



(2) 



which will be used frequently. 

For a given solution to the Bethe ansatz equations, the 
corresponding energy eigenvalue of the Hamiltonian is 
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E = - K [r ] - 2 l[(l 
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Note that this expression is independent of the spectral 
parameter u which can be chosen arbitrarily. The for- 
mula simplifies considerably with the choice u = oj, by 



employing (g), which yields a polynomial form. How- 
ever, for the purpose of an asymptotic analysis in the 
Rabi regime, it is more convenient to choose u = 0, while 
for the Fock regime we use u = rj 2 . 

Asymptotics. We start our analysis with the Rabi 
regime where n 2 N « 1. From the BAE it is clear that 
r f" v a — > f as 77 — > 0, so that v a f=a ztn^ 1 . However, when 
77 = we know that the Hamiltonian is diagonalizable by 
using the Bogoliubov transformation, from which we can 
deduce that the solution of the BAE corresponding to 
the ground state must have v a 1=3 rf~ . Excitations corre- 
spond to changing the signs of the leading terms in the 
Bethe ansatz roots. To study the asymptotic behaviour 
for the mth excited state, we set 



v a K,-ri l +e a +ri& a , a = l,...,m, 

v a ~ V _1 + e a + rj8 a , a = m + l,...,N, (4) 



with the convention that the ground state corresponds to 
m = 0. 

From the leading terms of the BAE for v a , a < m we 
find 



e Q - ep 



(5) 



which implies 



m(m — I) 



In a similar fashion we have for m < a < N 



N 



e a = - E Z — ' 



(6) 



which implies 

N N 

^ e Q = 0, E e " 

ot—m-\-l a— m+1 



(N - m)(N - to - I) 



It is clear from (g) and (|6|) why the Pauli exclusion 
principle applies in the present case. In the asymp- 
totic expansion for v a , e a is assumed finite. However, 
if v a = vp for some a, (3, then e a = ep and (||) and (j^) 
imply that e Q , ep are infinite which is a contradiction. 
Hence v a must be distinct for different a. Note also that 
for this approximation to be valid we require rj^ 1 >> e a . 
However, we see that \e a \ is of the order of N 1/2 . Thus 
our approximation will be valid for 77A 1 / 2 << I, which 
is precisely the criterion for the Rabi region and conse- 
quently N cannot be arbitrarily large for fixed 77, or vice 
versa. 

Now we go to the next order. From (0) we find 
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E_ m(m — 1) 
S a = 

a=l 



m(m — N) muj 2 



N 



E 5° 

a— m+1 



(N - to) (AT - m - 1) m(m-N) 



(N - m)ui 2 



which using (||) leads us to the result 

E m V 2 u; 2 (N-2m) V 2 N rj 2 

py — A+2m to A*— to). 

k 2 4 2 V 7 

The energy level spacings A m = £^ ni — S m _i are thus 



A m w k ( 2 + rfuj 2 + — (AT - 2m + 1) 



One may check that A m /N is of the order of A^ -1 . This 
indicates that the Rabi regime is semiclassical ||. Our 
value for the gap between the ground and first excited 
state agrees, to leading order in i] 2 N, with the Gross- 
Pitaevskii mean-field theory |ll| which gives a Josephson 
plasma frequency of uj = 2k(1 + ?] 2 N/2) 1 / 2 . 

Now we look at the asymptotic behaviour of the BAE 
in the Fock regime rj 2 >> N. It is necessary to distin- 
guish the following cases: (i) u = and (ii) w^O. 

(i) ui = 0. In this case, it is appropriate to consider 
the permutation operator P which interchanges the la- 
bels 1 and 2 in (|l|). For ui — 0, P commutes with the 
Hamiltonian, and any eigenvector of the Hamiltonian is 
also an eigenvector of P with eigenvalue ±1. Therefore 
the Hilbert space splits into the direct sum of two sub- 
spaces corresponding to the symmetric and antisymmet- 
ric eigenfunctions. From now on we restrict ourselves 
to the case when A~ is even, i.e., A" = 2M, although a 
similar calculation is also applicable to the case when A" 
is odd. A careful analysis leads us to conclude that the 
ground state lies in the symmetric subspace. The asymp- 
totic form of the roots of the BAEs for the ground state 
takes the "string" -like structure 



v j± 



-3> 



"(2j-l) 



+ M(M + l)rT 



j = !,••-, M. 



where C 3 M is a binomial coefficient. For this asymptotic 
ansatz to be valid, we require that any term in the asymp- 
totic expansion should be much smaller than those pro- 
ceeding. This yields i] 2 >> N which coincides with the 
defining condition for the Fock region. Throughout, the 
Pauli exclusion principle has been taken into account to 
exclude any possible spurious solutions of the BAE. 

The above structure clearly indicates that in the 
ground state the A" bosons fuse into M "bound" states 



and excitations correspond to a breakdown of these 
bound states. Specifically, the first and second excited 
states correspond to the breakdown of the bound state 
at — (M — l)rj, with the first excited state in the anti- 
symmetric subspace and the second excited state in the 
symmetric subspace. Explicitly, we can write down the 
spectral parameter configurations for the first two excited 
states 



vi+ w -Mi] + a 1+ r) 3 , vx- « — (M — 1)77 + ai-t] 3 , 
v j± w — (M - j)n + a j± r,-W-V , j = 2, • • • , M, 

with 



a 1+ 



02± 



M + 1 



<2i 



M(M + 1) 



-(M - l) 2 ± (M - 1)V13M 2 + 1QM + 1 
12 : 



a 3± w ± 



(M - 1)(M - 2) V / 2M(M + 1) 



24 

M-j + 1 



=%-i,±, J 



3,---,M, 



V(j + i)jC/-i)0'-2) 

for the (antisymmetric) first excited state and 

(M + 1)(2M+1) M(M + l) 

a l+ ~ n 5 a l- ~ 7, 



a-2± 



■(M - l) 2 ± i(M - 1)^/IIM 2 + UM - 1 



12 



03± ~ ±i 



a J± 



(M - 1)(M - 2) V / 2M(M +1) 



24 



M-j + 1 



^{3 + mj-m-2) 



aj-i,±, j = 3,--,M, 



for the (symmetric) second excited state. The breakdown 
of the bound state at — (M — j)r), j = 2, • • • , M results in 
the higher excited states. 

Substituting these results into (||) leads us to the 
asymptotic ground state energy 

E a fa -2Krf 2 M(M + 1), 

while for the first and second excited states we have 

, M 2 + M -2 



Ei fa nrj 2 — nrj 2 



2 o 5A/ 2 + 5M + 2 
E2 ~ Kf] + nrj . 



In contrast to the Rabi regime, the Fock regime is not 
semiclassical, as the ratio of the gap A and N is of finite 
order when N is large. 

We can perform a similar analysis for odd N . In this 
case, the gap between the ground and the first excited 
states is proportional to ni]^ 2 instead of nrj 2 . This indi- 
cates there is a strong parity effect in the Fock regime. 
Its possible physical implication remains to be explored. 
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(ii) w^O. In this case the root structure is somewhat 
more complicated than for u = 0, so we will not present 
the details. We remark however that our calculations 
show that up to order r\~ 2 the ground state energy eigen- 
value takes the same form as in the case u> — 0. Actually, 
the leading contribution arising from the ui term appears 
only as uj 2 r/~ 4 . This means that the results presented 
below are applicable for all values of uj (or equivalcntly 

a m ). 

Quantum fluctuations and coherence factor. Although 
it is difficult to define rigorously ||, the relative phase 
between Bose-Einstein condensates has been useful in un- 
derstanding interference experiments [|[i|]7|]24|l . By defini- 
tion, the relative phase <& is conjugate to the relative 
number of atoms in the two condensates n = Ni — N%. 
Using the Feynmann-Hcllman theorem, we find that 



< An 2 



>= 



dE 
OK 



4( 



dE 



For the ground state in the limit of strong tunneling 
(i.e., Rabi regime), < An 2 >= N - (A^N/Ej) 2 . In the 
case of weak tunneling (i.e., Fock regime), < An 2 >rj 
2N(N + 2){£j/K) 2 . The degree of coherence between 
the two BEC's can be discussed in terms of [|| 



1 

2N 



< a[a 2 



>= - 



1 dE 



where the last identity follows from the Feynmann- 
Hellmann theorem. In terms of the relative phase 
a =< cos$ > In the strong coupling limit, 

awl — N^ 1 (A^l) 2 /(8£j) 2 , indicating very close to full 
coherence in the ground state. In the opposite limit, we 
have a » 2(N + 2)£j/K << 1, indicating the absence of 
coherence. 

The above results give the first order corrections to the 
results presented in |6pl|] for the number fluctuations and 
the coherence factor at zero temperature. Moreover, we 
can extend this approach to gain insight into the effects of 
thermal fluctuations. This problem was recently studied 
in using Gross-Pitaevskii theory. In the Fock regime, 
the energy gap is large, and so for the low temperature 
thermodynamics it is sufficient to consider the first two 
excited states, for which the calculations are straightfor- 
ward. However, in the Rabi regime the small energy scale 
requires us to consider the entire spectrum. The asymp- 
totic form of the partition function is found to be (valid 
for all temperatures) 



m 

sinhf(iV + l) (3K 



sinh ^ 



4£j 



exp(/3N)T> 



l-exp/?(A + l) 
1 — exp 



with (5 = — (we set Boltzmann's constant ks = 1) and 
T> is the differential operator 



V 



(An) 2 N 
i£.jK 



N 

T 



V 1£jK 



N 
~2 



d_ 

8(3 



1 d 2 



In the usual way it is possible to derive a number of 
physical quantities from Z . For example, the tempera- 
ture dependent coherence factor is simply 



a(T) 



T d\nZ 
N d£j 



and an expression for the particle number fluctuations 
can be found in a similar manner. This provides an ana- 
lytic means to investigate the thermal effects for the Rabi 
regime. In order to investigate the crossover to weaker 
tunneling, it is possible to solve the BAE numerically 
to compute these quantities. This will be the subject of 
future work. 

In conclusion, we have shown that the canonical Hamil- 
tonian describing quantum tunneling between a pair of 
BEC's is integrable in the sense of the algebraic Bethe 
ansatz. An closed form expression was given for all of 
the energy eigenstates and eigenvalues in terms of the 
solutions of the Bethe ansatz equations. Our approach 
makes the calculation of form factors |20| and correla- 
tion functions possible without any of the limitations as- 
sociated with Gross-Pitaevskii mean-field theory or the 
number-phase representation. 
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